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ABSTRACT 

The 'shock in jet' model for cm-waveband blazar variabihty is revisited, ahowing for arbitrary shock 
orientation with respect to the jet flow direction, and both random and ordered magnetic field. It 
is shown that oblique shocks can explain events with swings in polarization position angle much less 
than the 90° associated with transverse structures, while retaining the general characteristics of out- 
bursts, including spectral behavior and level of peak percentage polarization. Models dominated by 
a force-free, minimum energy magnetic field configuration (essentially helical) display a shallow rise 
in percentage polarization and frequency dependent swing in polarization position angle not in agree- 
ment with the results of single-dish monitoring observations, implying that the field is predominantly 
random in the quiescent state. Outbursts well-explained by the 'shock in jet' model are present during 
7-ray flaring in several sources, supporting the idea that shock events are responsible for activity from 
the radio to 7-ray bands. 

Subject headings: galaxies: jets — magnetic fields — polarization — radiation mechanisms: nonther- 
mal — shock waves 



1. INTRODUCTION 

One of the early results from EGRET was the discovery 
that a subset of very active extragalactic objects, blazars, 
emit at 7-ray energies. A close association between ac- 
tivity in the radio band identified using total flux den- 
sity monitoring observations and detections by EGRET 
in the GeV band was e stablished in the mid-1990s, e.g., 
iValtaoia fc Terasrantal (1995), an association which has 
been confirmed using the large data base of measure- 
ments provided by Fermi and monitoring measurements 
from both single dish and V LBA imaging m e asure ments, 
e.g., iKovalev et aD ()2009( ): iRichards et all (|2010l ). For 
the past 25 years, the accepted explanation for blazar 
variability in the optical-to-radio bands has been the 
'shock in jet' model; hence the temporal associations 
between the activity in the radio and 7-ray bands and 
the fact that the EGRET detections occurred during 
the rise portions of the radio flares suggested to inves- 
tigators that the same shocks producing the radio flares 
were responsible for the emission in the 7-ray band, e.g., 
IValtaoia fc Terasrantal (1996). This shock scenario has 
been echoe d in recent s tudies of Fermi observations of 
AGN, e.g., lAbdo et all ()2010b[ ). but rigorous tests have 
not been extensively carried out. 

Evidence in support of the 'shock in jet' scenario orig- 
inally came from model fits to the broadband spec- 
tral evolution in the o ptical-to-radio band in 3C 273 
([Marscher fc Gea"illl985[) and independently from radia- 
tive transfer model fits to multifrequency total flux den- 
sity and linear polariza tion monitoring measuremen ts 
in the centimeter band (|Hughes. Aller fc AU"ei1ll989ai rbl. 
[T99I . The generally accepted scenario is that outbursts 
result from instabilities which develop naturally within 
jet flows, producing shocks. The magnetic fleld is ini- 
tially random within this emitting region pones et al.l 
I1985D . but the shocks produce a compression and an in- 
creased ordered component; the expected signature in 
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the radio band linear polarization light curves for such 
a shock event is a swing in the electric vector position 
angle (EVPA; an orientation orthogonal to the magnetic 
field direction in a transparent source) and an increase 
in the fractional linear polarization. This early model- 
ing assumed that the shocks had a specific orientation 
" transverse to the flow direction. Attempts to fit later 
radioband events in these same sources with the same 
model parameters used successfully in the original fits, 
however, failed, and the characteristic behavior of the 
variations suggested that shocks are more generally ori- 
ented obliquely to the jet flow direction; thus, the ini- 
tial transverse shocks identifled from their well-resolved, 
distinct appearance represented a special case of a more 
ge neral phenomenon. Conical shocks have been discussed 
bv lCawthorne fc CobbI (|1990l ). and the formulation pre- 
sented there in has been applied to the data for some 
sources, e.g., ' Marscher fc JorstadI ()2010[ ). 

Support for the shock model of major outbursts seen in 
single-dish data, and the propagating components seen 
in maps of parsec-scale flows, plus support for a shock 
explanation for at least some Fermi events, would come 
from a) 'revalidating' the 'shock in jet' model, by show- 
ing that oblique shocks can indeed explain the commonly 
observed reduced swing in EVPA through only tens of 
degrees, and associated increases in both fractional po- 
larization and total flux density (flares) with the spectral 
behavior exhibited in the data; b) showing that 7-ray 
flares occur at the same time as radioband events plausi- 
bly explained by an oblique shock model. It is also desir- 
able to test whether the single-dish data can discriminate 
between models involving purely random magnetic field, 
purely ordered magnetic field, or a mixture of the two. To 
those ends it is illustrated here that the salient features 
identified in radio band linear polarization and total flux 
density observations obtained by the University of Michi- 
gan monitoring program (hereafter UMRAO) are repro- 
duced by simulations using plausible input parameters in 
modeling involving radiative transfer through flows with 
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propagating, oblique shocks; and it is demonstrated that 
this shock signature is present during 7-ray flaring in 
several sources. 

In f|2]the results from the Michigan variability program 
are summarized, the general characteristics of the cm- 
waveband variability are defined, and their relation to 7- 
ray flaring is discussed, with examples. In ^the details 
of an oblique shock model for this cm- waveband activity 
are set out: the quiescent flow, models for both random 
and ordered magnetic field components, the structure of 
an oblique shock, and its propagation. The essential el- 
ements of an existing code for the transfer of polarized 
radiation are recapped in 21 a-nd the results of applying 
this to a fiow with representative flow parameters are 
presented in fJSl results for a range of observer orienta- 
tions and magnetic field topologies are contrasted. Sec- 
tion[6]summarizes and discusses the findings that oblique 
shocks can explain the observed cm- waveband events, the 
behavior of which support a model with at most a modest 
contribution from an ordered magnetic field component. 

2. OBSERVATIONS 

2.1. The UMRAO Database 

As part of the Michigan variability program carried 
out with the 26-meter radio telescope, source-integrated, 
single-dish total flux density and linear polarization mea- 
surements have been obtained for hundreds of sources, 
extending over time windows of up to 4 decades. The 
data acquired provide a comprehensive view of the prop- 
erties of blazar variability in the centimeter band based 
on long term data acquired and reduced in a consistent 
and uniform manner. These measurements commenced 
at 8 GHz in the mid-1960s, at 14.5 GHz in 1974, and 
at 4.8 GHz in 1978. The number of sources included 
in the UMRAO monitoring program increased dramat- 
ically in the late 1970s when an automated observing 
program was implemented; subsequently the fraction of 
time utilized per year for AGN observations increased 
to over 90%. Typically 20-25 sources are observed in a 
24 hour period; each observation consists of a series of 
on-off measurements over approximately 30-45 minutes. 
Measurements of calibrators are interspersed with ob- 
servations of program sources every 1 1/2 to 2 hours. In 
general, the observing cadence adopted for each program 
source is adjusted to follow the variability in that partic- 
ular AGN. Very active sources such as those included in 
the shock program described here are typically observed 
1-2 times per week at 14.5 GHz and once per week at 8 
and 4.8 GHz during flaring periods. This cadence is re- 
quired to follow the rapid variations in linear polarization 
expected from past observations. 

2.2. Characteristics Of The Variability 

Studies of the data for both individual sources and for 
statistical samples have identified a number of general 
properties of the variability: the characteristic variabil- 
ity time-scales of the total flux density, identified from 
first-order structure analyses, are of the order of two 
years in the observer ' s fram e with some spread, e.g., 
iHughes. Aller fc Alleii (|1992D . The ratios of peak to qui- 
escent amplitude of these variations are characteristically 
only a factor of a few, reaching at most a factor of 8. 
While numerous attempts have been made to identify 



periodic behavior in the total flux density light curves, 
in general consecutive events do not repeat. The out- 
burst shapes are characterized by a rapid rise followed 
by a more gradual decline. In the millimeter band, the 
Metsahovi monitoring data have been successfully fit 
assuming a function with an exponential rise, a sharp 
turnover, and an exponential decay with a time-scale 
approximately 1.3 tim es longer than the rise time-scale 
(jValtaoia et al.lll999D . but a generic shape is not appar- 
ent in the events observed in the lower frequency UM- 
RAO data, possibly due to the blending of individual 
events. The total flux density spectra in the centimeter 
band are typically inverted and relatively flat {\a\ < 0.5); 
however, these differ from source to source and can differ 
from event to event in the same source. The spectrum 
itself in many cases evolves during an outburst; in many 
large events the spectrum is inverted until burst maxi- 
mum, and flattens during outburst decline. Associated 
flares are often apparent in polarized flux density; and 
in some cases individual events can be seen in linear po- 
larization which appear as blended events in total flux 
density. 

Generic properties of the variability in total flux 
density and linear polarization are summarized in Ta- 
ble [TfA) , for comparison with models to be presented in 
Sjni The values listed are based on examination of the 
UMRAO database, while Figures [T][3] present specific ex- 
amples illustrating the range in these properties during 
recent well-resolved events temporally associated with 7- 
ray fiares detected by Fermi. The table contains: in col- 
umn 2, a simple measure of the flare shape given by the 
rise time, i.e. the time from start to peak, Trise, divided 
by the event duration, T; in column 3 the spectral in- 
dex at outburst start, at 14.5 GHz peak, and at end, in 
the UMRAO frequency range 4.8 to 14.5 GHz; in column 
4, characteristic values of the maximum fractional linear 
polarization (note that while the values listed are typical, 
values near 18% have occasionally been observed during 
fiares) ; in column 5 the monotonic swing in EVPA during 
an event, AEVPAtj„ie; and in column 6, the spectral vari- 
ation in EVPA throughout this swing, AEVPA/^eg- The 
spectral indices follow the sign convention S^, oc W^" , and 
the values listed for this parameter are based on numer- 
ical results for selected events, using paired values of 2- 
week averages of the data at 14.5 and 4.8 GHz. For linear 
polarization (peak fractional linear polarization and vari- 
ation in EVPA with time or frequency) , the values listed 
were determined from visual inspection of the long term 
UMRAO data. No corrections were applied for Faraday 
rotation; but events in sources with low source-integrated 
Faraday rotation measure were preferentially examined. 

Figure [T] shows the UMRAO monitoring data for the 
QSO PKS 1510-089. A series of 7-ray fiares occurred 
during the t ime period shown w hich hav e been analyzed 
in det ail in lAbdo et"ari (|2010a[ ) and in iMarscher et al.l 
(|2010( ). A particularly large and isolated event was iden- 
tified with st art date 10 March , 2009 and stop date 9 
Aprfl, 2009 (|Abdo et al.l l2010aD . The UMRAO light 
curves reveal that a radio band total flux density flare 
commenced in early 2009 at which time both a resolved 
linear polarization flare and an ordered swing in EVPA 
are apparent at 14.5 and 8.0 GHz. Because of the dif- 
ferences in variability time-scales in the radio and 7-ray 
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bands and the expected time delays produced by self- 
absorption in the emitting region, there are some uncer- 
tainties in associating specific events across bands; an 
event with the signature of an oblique shock, however, 
unambiguously occurred during this time period which 
includes the 7-ray flare. 

Figure [5] shows the radio band data during the period 
2008.0 through 2009.5 for the very active BL Lac ob- 
ject, OT 081 (1749+096), one of the three sources orig- 
inally modeled using the tra nsverse shock formulation 
(jHughes. Aller fc Ailed Il991| ) . Several resolved events 
are apparent in both the total flux density and frac- 
tional linear polarization light curves; an ordered swing 
in EVPA at 14.5 GHz occurred from August to December 
2008 followed by a shorter-duration ordered swing dur- 
ing January- February 2009. The Fermi light curve for 
this source is included in the variability study ba sed on 
data fr om the first year of Fermi operation ( A bdo et al.l 
l2010bD . A very large flare in the 7-ray band was in 
progress in August 2008, and a second flare which peaked 
in mid-March 2009 is also apparent in these data. Activ- 
ity, including strong flaring, occurred in both the 7-ray 
and radio wavebands. 

Figure |3] shows both the 7-ray and radio band light 
curves for O J 287. In this source the signature of a shock 
temporally associated with the outburst observed by 
Fermi is clearly apparent. The swing in EVPA through a 
limited range of about 40° is consistent with the passage 
of an oblique shock. 

3. JET MODEL 
3.1. The Quiescent Flow 

All computations are performed in a box 61 x 61 x 600 
cells in extent, with the axis of the jet parallel to the long 
axis, z, and polar angles are measured in the usual sense: 
9 from the z-direction, and (j) from the x-direction. A 
simple conical form has been taken for the jet boundary 
as observations cannot constrain more complex profiles 
explored over at most tens of jet radii. The opening angle 
(/z) is determined by setting the jet radius to be Ho on 
the inflow plane z = 0, and rhi on the outflow plane 
z = Zniax = 600; angles are typically 1 — 2°. 

Quiescent flow values are established across the entire 
domain, and then an apodizing filter 



a (p) = 0.5 



tanh 



-p + r (z) 



(1) 



1/2 

is applied, where p = (x^ + j/^) is the radial location 
of a cell at some z, where the jet radius is r (z), and ay, 
is a parameter that controls the flow's boundary extent; 
fltu is typically one cell. 

Whether flows exhibit significant acceleration or de- 
celeration remains controversial, and indeed, might dif- 
fer from source to source. D eceleration has been in- 
ferred for FR 1 radio galaxies (jLainei [19961 iLaing et"al] 
Il999( ). while a recent study of the BL Lac object OJ 
287 finds that the flow remains highly relativistic to dis- 
tances as high as hund reds of kiloparsecs from the nucleus 
(|Marscher fc Jorstadl[201Q) . In view of this uncertainty, 
and the limited number of jet radii explored, the quies- 
cent flow speed is taken as a constant specified by its 
Lorentz factor, 7/; values are typically 2 — 5. A diverging 



flow is then established, with stream lines parallel to the 
z-axis on axis, and parallel to the flow boundary there. 
The magnetic field is constructed as described in 
and as s umed to fall in strength along the flow as r (z)~^. 
iLeahvl ([1991 gives a good overview of the evolution of 
parameters in a diverging adiabatic flow, and it might 

be thought that a fall-off as r (z)^^^^ would be appro- 
priate for a flow of fixed speed, given that the default 
assumption is that of a random field. However, that as- 
sumes no coupling between the components, so that the 
component perpendicular to the fiow falls more slowly, 
driving the flow away from isotropic turbulence. The 
assumption of an isotropic random component at every 
location in the quiescent flow implies a turbulent driving 
that persists the length of the flow, transforming per- 
pendicular into parallel field, which declines with the as- 
sumed r (z) ^ dependence, and vice versa, coupling the 
decline of the perpendicular component to that of the 
parallel component. Further, the helical ordered field 
that is also explored, if generated from a turbulent com- 
ponent by a dissipative dynamo process such as noted 
in ! j3.2.2l will then follow the same global trend. The 
adopted dependence is thus a simplification, but reason- 
able given the other approximations. As discussed in the 
next section, the method used to generate a random field 
component, the ad hoc scaling along the fiow, and the 
scaling of the relative strengths of random and ordered 
field components, mean that the field does not satisfy a 
self-consistent magnetohydrodynamic model for its gen- 
eration, and is not divergence-free. However, the aim 
is to construct a plausible magnetic field topology with 
which to explore the properties of the emergent radiation, 
not a field to be evolved subject to constraints to ensure 
conservation and/or preservation of the divergence-free 
character of the field. 

The fiow is assumed to be pervaded by a power law 
distribution of radiating particles, with density 



n (7) d-y = no7 ^ d-y, 7 > 7^, 



(2) 



with the index S fixed. Adiabatic gains and losses will 
not change the slope, synchrotron losses are ignored over 
the tens of jet radii propagation explored here, and no al- 
lowance is made for Fermi acceleration at shocks, whose 
primary role will be to produce a small number of par- 
ticles radiating at frequencies above those probing the 
quiescent emission. The absolute scaling of the density 
is arbitrary, as functions of the density constant no, and 
the magnetic field strength are absorbed into a fiducial 
optical depth, r, which acts as a free parameter; see [jH 
The density constant, uq, is taken to fall due to adiabatic 
expansion of the fiow according to r (2;)"^^'''''^''^'^ , and the 
cutoff thermal Lorentz factor (see fJ3| for the significance 

of this quantity), 7;, falls as r (z) "^^^ . 

3.2. Magnetic Field Structure 

Arguments in favor of a largely tangled magnetic field 
were made by iHugh esi (2005), the salient points being: 
a) the low degree of polarization exhibited by compact 
extragalactic radio sources when in a quiescent state has 
been widely interpreted as due to 'root-N' depolariza- 
tion in a synchrotron source with many randomly ori- 
ented magnetic cells; b) there have been successful mod- 
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els of the temporal, spatial, and spectral attributes of 
outbursts in a number of individual sources, with a sce- 
nario in which the shock compression of a flow provides 
an effective order to the otherwise random magnetic field; 
c) the observed levels of linear and circular polarization 
are best modeled by a scenario in which root-N depolar- 
ization plays a significant role, albeit that some of the 
magnetic energy is in an ordered component. 

In contrast to this pictu r e, Gabuzda and coworkers 
jGabuzda et al."2004', '2005*; Mahmud & Gabuzda '200^; 
lO'SuUivan & Gabuzda 2008; Contopoulus. ct al. 2009) 
have argued in favor of a large-scale magnetic field - in 
particular with a helical character - the key evidence 
being rotation-measure gradients across jets, which are 
interpreted as due to the systematic change in the line- 
of-sight component of the jet magnetic field. Further- 
more, it is almost invariably the case that 'central en- 
gine' models invoke an ordered magnetic field in the en- 
vironme nt of a supermassive bla c k hole (see, by way of 
exam ple iKomissarov et afl 120071 : iMcKinnev fc NaravanI 
120071 ) and this order (together with that due to shear) 
would be expected to be imprinted upon the parsec-scale 
flow. Recent work by Brodcrick & McKinncy (2010) and 
iTavlor fc Zavalal (|2010[ ) notes that the presence of a he- 
lical field does not necessarily imply that a simple mono- 
tonic variation in rotation measure across the flow will be 
observed, and the resolution currently available precludes 
definitively establishing the character of any ordered field 
within the jet. However, clearly there is evidence consis- 
tent with such a field in a number of sources. 

If an ordered field exists, and its sense is determined 
by the spin of the central black hole, and/or the rotation 
of the associated accretion disk, then the sign of the re- 
sultant circular polarization would be expected to be a 
fixed attribute of any source. However, Stokes V mon- 
itoring does not provide a simple picture; while, some 
sources exhibit such a constancy, in others changes in 
the handedness of circular polarization with time and fre- 
quency are qui te unambiguous ov er durations of months 
to years, e.g., lAUer. Aller. fc Hu ghes (2010). Also, ev- 
idence for a reversal in the direction of the rotation- 
m easure gradient across the jet of 1803-1-784 presented 
bv lMahmud et al.l ()2009() . at large distance from the core 
and away from the portion of the jet most affected by fi- 
nit e beam size, are at odds with t he results of simulations 
bv lBroderick fc McKinnevI (|2010( ). who argue against the 
'magnetic tower' model for a temporal change in the sign 
of the rotation measure. The analysis of all-S tokes po- 
larization data for 3C 279 (|Homan et al.|[2009l) also does 
not find support for the role of a helical field in the 
production of circularly polarized emission. Some res- 
olution of these conflicting perspectives may yet come 
from a more careful consideration of the emission loca- 
tion (the jet spine) and the source of Faraday rotation 
(the jet sheath; the do minant contribution to Farad ay 
effects in the models of iBroderick fc McKinnevI (|2010[ )). 
Until higher spatial resolution data become available, 
and simulations can track the distribution of radiating 
particles, the topology of the magnetic field remains un- 
certain, and a contentious issue, and so here both random 
and ordered fields are considered. 

3.2.1. Random Field Component 



An elegant method of establishing a random magnetic 
field is to select random phases, and random amplitudes 
from a Rayleigh distribution, for the Fourier transform 
of the vector potential, A(k). The Fourier transform 
of the magnetic field is then B (k) = zk x A (k) , and 
an inverse Fourier transform of this yields a magnetic 
field that is guaran teed to s a tisfy t he divergence- free con- 
straint. (See, e.g.. iTribbld (|1991l ) for a detailed exposi- 
tion of thi s approach.) However, as first noted in this 
context bv lJonesI (|1988[ ). such a technique leads to large- 
scale Fourier components that give rise to a small but 
significant excess polarized flux density. An exploration 
of this in the current study showed that unacceptable 
levels of polarized flux density were produced, and that 
t here is no s imple way to mitigate that effect. 

iJonesI ()1988i) used a lattice of nested cells, each with 
randomly chosen field orientation, and with the ampli- 
tude of the fleld components scaled to produce an ap- 
proximately Kolmogorov power spectrum. The goal of 
the current study is such that there is a need to exten- 
sively explore a very large parameter space, and thus use 
a 3D grid of limited resolution, in order to keep individ- 
ual runs to a modest time. Furthermore, a jet length 
of many jet radii is required, because although the total 
flux density falls rapidly after a propagating enhance- 
ment to the particle and field densities passes the r = 1 
surface, a polarized fiux density significantly above the 
quiescent value persists for much longer. Thus there are 
too few cells within the jet to effectively employ a multi- 
scale approach, and simply assign randomly chosen field 
directions within each cell. In principle the field should 
be advected with the fiow, but such advection cannot be 
handled self-consistently at oblique shocks, due to the 
simple way that those structures are modeled (see §3.4|) , 
and as this study is concerned with 'macroscopic' behav- 
ior in the polarized emission, not random fiuctuations 
arising from fiow structures lighting up different realiza- 
tions of random field, a static random field is built into 
the jet. 

To reduce the computational overhead, multiple real- 
izations of a random field have not been generated for 
each parameter combination. However, in ^15. 21 radia- 
tion transfer is performed for a range of azimuthal an- 
gles about a flow with propagating oblique shock. By 
symmetry the results for azimuthal angles 225°, 270°, 
and 315° are identical to results for 135°, 90°, and 45° 
respectively, except that the different angles sample dif- 
ferent projections of the random fleld component, and so 
the differences in emergent radiation for the three pairs of 
runs provide insight into variations to be expected from 
different realizations. At all azimuthal angles, total in- 
tensity curves are identical, as are those for percentage 
polarization, but during outburst there are small differ- 
ences in EVPA as a function of frequency; this is greatest 
at the first time shown, the quiescent state, because of 
the low level of polarized emission then. 

3.2.2. Ordered Field Component 

As discussed above, ordered, and in particular 'heli- 
cal', magnetic fields have received much attention over 
the last decade. In order to include the most physically 
plausible ordered field component, with the least number 
of free parameters, a force-free, minimum energy config- 
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uration is adopted. On the sub-parsec scale of interest 
here, the flow is far enough from the central engine that 
local evolutionary effects should dominate over the influ- 
ence of the engine's ergosphere or inner accretion disk. 

Indeed, long-term changes in the character of a large- 
scale magnetic field are possible, if jets are subjected to 
turbulence and a dissipative dynamo process. Turbu- 
lence could be driven throughout the body of the flow 
due to their ultra-high Reynolds number. Mean field dy- 
namos in s heared, turbulent flows are known to occur 
(jRogachevs kii & Klccorin 2003), and although the pro- 
cess has not been explored in cylindrical geometry, it 
can be speculated that the mean field then takes a form 
similar to that of a force-free flux tube. This field ex- 
ist s in what would be the asymptotic region discussed 
bv lMcKinnevI 1)20061 ) . and its symmetry is thus unrelated 
to the angular momentum of the 'central engine'. Al- 
though the process of field reversal is different from that 
exhibited by the solar-terrestrial field - bei ng stochastic 
rather than periodic - recent simulations (jYousef et al.l 
|2007|) have shown such behavior, and it would naturally 
explain the observed circular polarization 'flips'. 

Whether the fleld evolves due to dynamo action or 
relaxation, a likely end point is a force-f ree configura- 
tion V X B = gB (lEilek fc H ugheillQQlD. It has been 
shown bv lKonigl fc Choudhu ri (1985) that in cylindrical 
geometry, only the modes m = and m = 1 contribute 
to the general solution to the force-free equation in the 
minimum energy state. For a range of jet and ambi- 
ent medium parameters, and for any parameters leading 
to a jet with small deviation from axisymmetry, only 
the 171 = mode need be considered. Then, by the 
divergence- free constraint V • B = 0, the axial wavenum- 
ber k = 0. 

In circular cylindrical coordinates, {p,(j),z), this leads 
to a field configuration 

B = So (0, Ji {Kp) , Jo {Kp)) , (3) 

where the J are Bessel functions of integer order, and K 
is a constant that must be determined. The only phys- 
ical boundary condition is that the radial component of 
B vanish there, but as this is identically zero for the 
m = mode, this boundary condition provides no con- 
straint. By choosing K = 2.41/r(z), with r {z) the jet 
radius, the axial fi eld goes through a null at the boundary 
(|Lundquistlll950( ) . While it is in principle permissible for 
the axial field to reverse sign inside the jet, the adopted 
value prevents that, and is consistent with the simple 
field topology selected by evolution to a force-free, mini- 
mum energy state. This field is adopted 'locally', in that 
the jet radius increases due to fiow divergence, but note 
that it is not built self-consistently into the prescription 
of the magnetic field. 

The constant Bq is varied to model fiow divergence and 
shock compression. The strength of the random com- 
ponent is established by scaling so that < B^^^^^ > is a 
specified multiple of the ordered field strength. In the 
absence of a fully self-consistent magnetohydrodynamic 
model, it is not clear whether the random field strength 
should be a multiple of the local ordered field strength 
(corresponding, for example, to a turbulent generation 
of the former from the latter by a fixed number of eddy 
turn-overs); or spatially fixed for a given jet radius, a 
multiple of the axial ordered field strength (correspond- 



ing, for example, to a turbulent generation of the former 
to fixed multiple of the kinetic energy density in fiow 
turbulence). Both approaches have been tried, revealing 
that the choice has no impact on the final results, merely 
changing the precise value of p = '^t-cdl < ^ran 1 
where a transition in the characteristic behavior of the 
polarized fiux density is evident. 

3.3. Axial Magnetic Field 

It is extraordinarily difficult to define a quiescent state 
from data. The UMRAO monitoring data exhibit almost 
continuous activity, and even apparently inactive phases 
may be periods in which small-amplitude outbursts are 
temporally unresolved. VLBI monitoring similarly re- 
veals few if any inactive epochs, and apparently-quiescent 
fiows may contain spatially unresolved components. The 
UMRAO database has been searched to identify periods 
of 'quiescence', during which there is measurable linearly 
polarized emission. Between 1998 and 2001 07354-178 
exhibited a low level of constant fiux density at the UM- 
RAO frequencies, with ~ 2% polarization and EVPA 
implying an ordered magnetic field lying within some 
tens of degr ees of the jet direc tion determined from MO- 
J AVE data ([Lister et al.ll2009l ) . Similar numbers apply to 
NRAO 530 from 2002 to 2010 as iUustrated in Figure H 
Based on (admittedly sparse) examples such as these, it 
is plausible to assume that in the quiescent state a weak 
axial field is added to any random component, and in all 
simulations that follow, an axial mean field with 2% the 
energy density of the random component is added. This 
establishes a well-defined EVPA in the quiescent state. 

3.4. Oblique Shocks 

Modeling a transverse structure comprising forward 
and reverse shocks, separated by a contact discontinu- 
ity, is straightforward, as no lateral flow is required: the 
shocked flow domain expands in the frame of the con- 
tact surface, and either the expansion may be ignored 
during the brief interval that the structure traverses the 
r = 1 surface, or the expansion can be modeled self- 
consistently from the jump conditions. Furthermore, ei- 
ther the forward or reverse shocked flow can be assumed 
to dominate emission, or emission from both regions can 
be included. 

A kinematic model for an oblique shock is difficult 
to construct, because by its very nature, an oblique 
shock deflects the flow. Indeed, they have lon g been 
invoked to model flow curvature in radio jets ()Smithl 
Il984f ). Deflection of the flow would appear to be in- 
consistent with the simple, static, conical flow boundary 
assume d here. However, c onsider BL Lac as an example 
source: iDenn. Mutel fc Marscheii (j2000f ) present exten- 
sive monitoring by VLBI which reveals knots with lin- 
ear polarization behavior strongly suggestive of oblique 
shocks. While these authors stress the curvature of the 
knot trajectories, and difference between the kinematics 
of the four components studied, an overlay of the tra- 
jectories (as defined by the data, rather than the helical 
fiow model) for components S7 and S9 exposes a remark- 
able similarity, while the trajectories of S9 and SIO differ 
significantly only at the last epoch of the later, and in 
that SIO may be followed closer to the core. While com- 
ponent S8 does follow a distinctly different, and more 
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curved, trajectory, the general impression is of a well- 
defined channel with modest curvature amplified in pro- 
jection, that can support the propagation of multiple 
ob lique structure s. Furthermore, studies such as that 
bv iHughesI ()2005l ) show that oblique shocks (which may 
be transient features) can form and propagate, filling the 
cross section of a flow, while that flow maintains a simple 
jet-like, rectilinear form. 

An oblique shock is thus taken to be a plane spanning 
the jet, of obliquity 77, measured with respect to the di- 
rection of the upstream flow, which is approximated as 
being in the sense of the jet axis, ignoring the small di- 
vergence of the stream lines. 77 = 90° corresponds to a 
transverse shock, and the orientation of the shock is spec- 
ified by the azimuthal direction of the shock normal, ip. 
The shocked fiow is a volume extending from this plane 
to a parallel plane, a distance w along the flow axis. The 
limit of the shocked fiow is fixed; it is on the 'core' side 
of the shock plane for a "forward" shock (moving faster 
than the underlying jet flow), and on the other side for 
a "reverse" shock (being overtaken by the jet flow). The 
shocked flow may be terminated by adiabatic expansion 
of the shocked material, radiative losses within the flow, 
or simply reflect the distance the downstream flow has 
extended since formation. Adiabatic expansion is not 
relevant in this context in the current model, as by con- 
struction, the flow is constrained to lie within the quies- 
cent jet boundary. Radiation losses could have imposed 
a frequency-depen dent length on the downstream fiow 
(jMarscher fc GeaHflQSS,) prior to evolution through the 
modeled domain, but as explained in gl the UMRAO ob- 
serving frequencies are close (being separated by ~ -s/S) , 
and the length determined by radiation losses goes as 
so little frequency-dependence of the structure is 
expected. For the sake of definitiveness, it is assumed 
that the length of shocked flow is simply the evolved ex- 
tent of the downstream flow, which changes little during 
the brief propagation through the modeled domain. 

The shocked flow is characterized by a com pression 
K < 1, and following iCawthorne &: CobbI ()1990D . the up- 
stream flow speed in the rest frame of the shock, f3^ is 
computed: 



5'2. 



{k^ sin^ 77 - cos^ 77 - 1 -I- [(k^ sin^ 77 - cos^ 77 - 1)^ 
4/t^(8sin^ rj + l) sin^ 77 - 4cos^ v]^}/ 
{2[K^(8sin^ 77 H- l)sin^77- cos^ 77]}. 



(4) 



From lLind fc BlandfordI ()1985D it follows that the down- 
stream flow speed (also in the shock frame) is given by 



13d 



(^1 ~ p'^ cos^ rj^ + dp'J" cos^ r] sin^ 7/ 



3/3,^ sin 77 

and the flow is deflected by angle ( given by 



tan 



tanC 



tan 77 (tan^ 77 + 1 2/3^2) 



(5) 



(6) 



The flow direction within the shocked region is modi- 
fied by the shock speed, /3s, is computed from the as- 



sumed quiescent fiow speed {(3f) and the upstream speed 
in the shock frame (/?'«), and the downstream flow speed 

(3^ is transformed into the observer's frame, (3d- Within 
the shocked region of the jet, the particle density is in- 
creased by K^^'^+^^Z'^, and the magnetic field components 
projected on the shock plane are increased by the com- 
pression factor. 

3.5. Time Evolution 

As noted in ^'6.1l evidence for the acceleration or decel- 
eration of fiows is ambiguous, and comes primarily from 
the observed deviations from simple rectilinear motion 
displayed by a number of superluminal components. It 
is difficult to distinguish between an acceleration or de- 
celeration of the bulk flow, and that of shock waves on 
the flow, in addition to which acceleration or deceleration 
might be in direction only , due to a (mod est) curvature 
of the jet (see for example IMarscheil l2006). 

In the absence of compelling evidence for shock ac- 
celeration or deceleration, the propagation of the shock 
is modeled by a simple constant displacement of the 
shocked region discussed in i 33.4( radiation transfer is per- 
formed for the quiescent fiow at the first time step, and 
then over the remaining — 1 ^ 10 time steps the cen- 
ter of the shocked region traverses the flow so that the 
shocked domain lies just within the inflow jet volume at 
the second time step, and would have left the jet volume 
after the last time step. The corresponding physical time 
will be determined by the shock speed {/3s), computed as 
described in ^ 33. 41 and the physical length of the fiow; 
as the latter is subject to an arbitrary scaling, the time 
evolution is arbitrarily scaled to match typical outburst 
periods seen in the UMRAO data. 

4. RADIATION TRANSFER 

The transfer of polari zed radiation i s performed using 
the method described bv lHughei (|2005f ). and is briefly re- 
capped here. It is based on a formulation for the transfer 
of polarized radiation through a diffuse plasma, allowing 
for emission, absorption, the birefringence effects of Fara- 
day rotation and mode conversion (which can produce 
modest levels of circular polarization), and relativistic 
aberration and boost, which has been described in detail 
by a number o f authors, and is compactly summarized 
byQonei (p88^. 

The observer lies at arbitrary polar angles {0,</>), de- 
fined in the conventional sense with respect to the Carte- 
sian system used to describe the kinematic model. An 
array of lines-of-sight is established, using a preset den- 
sity of lines along the longest axis of the projection of 
the computational volume on the plane of the sky, with 
a commensurate number orthogonal to that, to ensure 
equal resolution in the two directions. For the results pre- 
sented below, a resolution of 128 pixels along the longest 
axis was used. For each line-of-sight, the algorithm finds 
the most distant cell, and radiative transfer is performed, 
cell to cell, until the 'near side' of the volume is exited. 
Within each cell, an aberrated magnetic field direction is 
computed from the rest frame field, velocity, and observer 
location, and is used with transfer coefficients modified 
by the relevant Doppler boost. For a given epoch of 'ob- 
servation' cell values corresponding to the appropriate 
retarded time should be used. However, as this leads to 
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an enormous increase in the computational overhead, and 
can lead to simulated VLBI maps that are very difhcult 
to interpret, because of the problem of unambiguously 
associating map features with physical flow structures, 
the radiation transfer is initially performed without re- 
tarded time effects included, and the latter will be con- 
sidered only when critical to understanding particular 
features of light curves or maps. Simulations to justify 
this strategy are presented in ii5.ll 

Radiation transfer is performed at dimensionless ob- 
serving frequencies of 0.6, 1.0, and 1.8125 corresponding 
to the U MRAO values of 4.8, 8.0 and 14 .5 GHz. As de- 
scribed in iHughes. Aller fc Alleii (|1989aD a fiducial ther- 
mal Lorentz factor (7c; 'thermal' implies the Lorentz 
factor associated with the random motion of the emitting 
particles, as opposed to the bulk flow Lorentz factor) is 
adopted, which is the energy of those particles radiating 
at the central observing frequency in the average mag- 
netic field for the jet volume at the first time step of 
evolution. The value of this fiducial Lorentz factor is not 
in itself significant, but model fitting to the data con- 
strains how far below this the power law electron energy 
distribution extends (to 7^), because it is the presence 
of these lower energy, but still relativistic, particles that 
produce Faraday effects in the absence of a 'cold' parti- 
cle distribution; adopting a plausible fiducial value thus 
enables us to estimate the low energy cutoff to the power 
law distribution. 

Prior to the radiation transfer, the optical depth for 
each cell must be known; this is computed from the 
dimensionless line-of-sight length through the cell, cell 
magnetic field, and cell particle density. This is scaled 
with a single adjustable parameter chosen to produce 
some 'target' optical depth through the entire volume at 
the central observing frequency, in the average magnetic 
field and particle density for the jet volume at the first 
time step of evolution. An initial 'target' optical depth is 
chosen in light of the spectral characteristics of the data 
being modeled; the value is then adjusted to fine-tune 
the model fit. Given the actual observing frequencies, 
the choice of fiducial Lorentz factor implies a particu- 
lar magnetic field strength; given that, and the optical 
depth needed to reproduce the data, knowledge of the 
physical scale of the jet implies a value for the parti- 
cle density, or vice versa. However, as the concern here 
is only with learning what such model fitting can say 
about the topology and orientation of observer, shock, 
and magnetic field, the implied field and particle densi- 
ties are not explored. 

5. MODEL RESULTS 

5.1. Role Of Time Delay 

As discussed in 21 our primary goal is to 'revalidate' 
the 'shock in jet' model, demonstrating that oblique 
shocks retain the temporal and spectral characteris- 
tics of the total and polarized fiux density behavior 
of cm-band outbursts previously modeled with trans- 
verse structures, while accommodating more complex be- 
havior of the EVPA. As 'proof of concept' parameters 
are chos en similar to those used to m odel outbursts in 
BL Lac (IHughes. Aller fc Alleri[l989bh . using a bulk fiow 

^ This docs not imply a thermal distribution of particle energies. 



of Lorentz factor jf = 2.5 and compression k = 0.7, but 
with a shock oriented at 45° to the flow axis, with the 
shock deflection in the plane of the sky; see Run A in 
Tabled 

Assuming a forward shock (F), implies the shock plane 
moving over the underlying fiow at Lorentz factor 6.7, 
which when viewed at an angle to the flow axis of 10° 
would produce an apparent superluminal motion of the 
leading edge of /3app ^6.5. This is a significantly faster 
shock speed than adopted in modeling the 1980s out- 
bursts, but the latter adopted a reverse (R) shock model 
with the observer at a larger angle to the flow sense, 
guided by the low apparent speeds reported in the liter- 
ature at the time. The current choice of parameters is 
guided by the much higher apparent speeds reported in 
the more recent literature (Lister et al. 2009), bearing in 
mind that over decades of activity the jet fiow direction 
can change, either forward or reverse shocks might be 
evident at different epochs, and that shock events within 
a single source could exhibit a range of obliquities. As 
shown in it is possible for oblique shocks to give rise 
to even higher percentage polarization, and thus more 
prominent, events than a transverse shock with similar 
parameters, for fiows seen within tens of degrees of the 
line-of-sight. 

Of the other parameters, optically thin spectral index 
is reasonably well constrained by observation and has 
been fixed accordingly. Others, such as fiow opening an- 
gle and fiducial Lorentz factor, were picked on the basis 
of typical values discussed in the literature over many 
decades, to see if choice of such 'plausible' values admits 
model light curves with the same characteristics as those 
observed. Cutoff Lorentz factor, shock width and optical 
depth (see ij4]), are adjusted in an attempt to reproduce 
light curves looking most like the example UMRAO data. 
The 'order fraction' and orientations are subject to study 
in later sections. 

Figure [5] shows the evolution of flux density, percent- 
age polarization and EVPA for this model, without the 
inclusion of retarded time effects. As described in SjHthe 
calculations are done using dimensionless quantities; to 
'guide the eye', time and flux density have been arbitrar- 
ily scaled to values representative of those seen in the 
UMRAO data. The typical characteristics of UMRAO 
bursts described in i)2.2| - the fractional flux density in- 
crease AS/ < S >, spectral evolution through a partially 
optically thick phase, percentage polarization with opac- 
ity/Faraday effects evident at the lowest frequency, and 
swing in EVPA by tens of degrees - are all reproduced. 
Note that in earlier modeling it was found that to fit 
the spectral characteristics of the polarized fiux density, 
a fairly low cutoff thermal Lorentz factor 20) was 
needed for some sources. Here the general character- 
istics of UMRAO outbursts are well-reproduced with a 
value that means opacity effects dominate, with Faraday 
effects being only marginal: this implies that in general 
few low energy electrons are present in these sources. 

While use of retarded time in the modeling would be 
necessary for detailed fits to data, the general character- 
istics of the total and polarized flux density light curves 
and spectral behavior, even without using retarded time, 
reproduce the behavior exhibited by the data f ^2.2p . Ta- 
ble[ljB) presents the model values corresponding to those 
derived from the UMRAO database, and discussed in 
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^2.2^ A comparison of these is made in SjHl 

5.2. Orientation 

Figure m shows models Bl, B2, etc. - the same model 
as in Fig. [5] (see Table [2]), but for a range of observer 
orientation with respect to the shock plane. Recall that 
ip = 0°, so that the shock normal lies in the x-z plane 
of the Cartesian coordinate system. The observer ori- 
entations explored here (0°, 45°, 90°, 135°, 180°, 225°, 
270°, and 315°) correspond to starting with a view par- 
allel to the cc-axis, and then rotating around the jet so 
that by panel (E) the observer is 'behind' the shock. 
As one would expect, the total flux density light curves 
are minimally changed by a change in azimuthal orienta- 
tion. For an observer orientation within tens of degrees of 
0obs = 0°, orientation and aberration conspire to provide 
a more nearly 'face-on' view of the shocked flow, and the 
percentage polarization is small. Indeed, at (j)obs — 0° 
the polarized flux density from the shock region cancels 
the small orthogonally polarized flux density associated 
with the axial field, as evidenced by the varied behavior 
in EVPA in panel (A). However, the percentage polariza- 
tion is approaching 5% by orientations of 45 and 335°, 
and exceeds 10% at peak over a large range of angles. 
The important conclusion is that azimuthal orientation 
does not play a major role in the total and polarized flux 
density outburst light curves; special conditions do not 
need to be invoked for polarizations of this order to be 
seen, and (subject to flow speed and polar orientation) 
most oblique structures will give rise to significant levels 
of polarized fiux density. As the azimuthal angle does 
not play a significant role, a value of 90° is adopted in 
what follows. 

5.3. Oblique versus Transverse Shocks 

In this section the run of percentage polarization with 
jet inclination for the transverse case is established, as a 
measure by which to judge the behavior in the oblique 
case; it is then shown that adopting an oblique shock 
does not radically change the behavior of percentage po- 
larization - in fact leading to slightly higher values at 
some angles. This provides further evidence that the 
'shock in jet' model survives the introduction of oblique 
structures necessary to explain the temporal EVPA be- 
havior seen in the UMRAO data, and the evolution of 
features found in time sequences of VLBI maps. 

Figure [7] shows models CI, C2, etc., ~ the same model 
as in Fig. [5] (see Table [5]) , but contrasting transverse and 
oblique shocks for a range of observer orientation with 
respect to the flow axis. As orientation is changed, the 
free parameter t is adjusted to ensure a similar spectral 
behavior of the total flux density. In panels (A-C) the 
shock is transverse, so the azimuthal location of the ob- 
server plays no role in determining appearance, and the 
observer is viewing at angles 20°, 40°, and 60° respec- 
tively to the jet axis. Note that with increasing angle the 
decline in percentage p olariz ation is rather slow. From 
iHughes. Aller fc Ailed (|1985l ). a flow with compression 
K = 0.7, seen edge-on, and in the absence of opacity and 
Faraday effects, would be expected to exhibit polarized 
emission ~ 25%, dropping to ^ 8% at an angle of 50°. 
Given the modest Lorentz factor of the shocked flow in 
the observer frame, ~ 3, radiation from this angle outside 



the critical cone of the flow (namely, at 140° to the flow 
axis in the flow frame) would be seen by the observer 
viewing at 50° to the flow axis, an orientation spanned 
by panels (B) and (C) in the flgure. In the simple trans- 
verse case, for this level of compression, quite high levels 
of polarization will be seen well beyond those values of 
viewing angle usually adopted in blazar modeling. Ta- 
ble [IJC) presents the model values for Run CI, corre- 
sponding to those derived from the UMRAO database, 
and discussed in ^ 32. 21 The only large difference com- 
pared with the values presented in Table [HB) is in the 
jump in EVPA as expected, as well as in the spread in 
EVPA at outburst end. 

In panels (D-F) the original oblique shock is viewed at 
angles 30°, 50°, and 70° to the jet axis, for an azimuthal 
orientation of 90°. (For an azimuthal orientation of 0° 
the situation will be approximately - subject to flow de- 
flection and different aberration - as for the transverse 
case, modulo an angular offset.) It can be seen that a 
high level of polarized emission persists to quite large an- 
gle from the flow axis, a result related to the geometrical 
effect that as the shock becomes more nearly parallel to 
the flow, rotation in polar angle has no effect on the ap- 
pearance of a flow seen initially nearly edge-on. As with 
azimuthal orientation, it is concluded that special condi- 
tions do not need to be invoked, and that the observed 
levels of percentage polarization can be seen for oblique 
shocks for a signiflcant range of observer orientation with 
respect to the flow axis. 

5.4. Magnetic Field 

Figure [5] shows models Dl, D2, etc. - the same model 
as in Fig. [5] (see Table [2]), but for a range of 'order mul- 
tiple', a measure of the importance of an ordered com- 
ponent of magnetic held. The values explored here (0.1, 
0.3, 1.0, and 3.0) correspond to moving from an essen- 
tially random held to an almost totally ordered field. A 
fixed multiple a t each point within the flow is assumed; 
as discussed in ii3.2.2[ trials show that the details of how 
the relative strength of these two components is modeled 
does not have a signiflc ant influence on the results. 

As noted bv Lvutikov et al.l (|2005j ). low levels of per- 
centage polarization may result even if the emitting vol- 
ume contains wholly ordered magnetic field, and that is 
evident in the values seen here outside of shock events 
- significant compared with the random field case, but 
low relative to the maximum value for optically thin syn- 
chrotron radiation. Careful choice of field topology might 
reduce the value further, but as argued in ^3.2.21 the 
adopted one appears to be the most physically plausible. 

The nature of the outburst in percentage polarization 
differs markedly between the case of polarization induced 
by compression of a random field, as seen in panel (A) of 
Fig. m and the case of polarization induced by asymmet- 
ric compression of an ordered field, as seen in panel (D). 
In the latter case the amplitude of the rise is limited (by 
less than a factor ^ 1.5), with a strong frequency depen- 
dence unlike either that seen in panel (A), or observed 
outbursts; the frequency dependence arises because at 
the lowest frequency opacity causes polarized emission 
to be dominated by counter-directed magnetic fleld arcs 
near to the jet boundary. The behavior of the EVPA 
is also markedly different between the cases. For the 
random field example the EVPA swings abruptly, and 
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similarly at all three frequencies, to an approximately 
constant value during the long decline in outburst. In 
the ordered field example there is a strongly frequency- 
dependent decline, followed by an immediate recovery, 
and no distinct plateau of EVPA during outburst. The 
evolution shown for the cases of a predominantly ordered 
(helical) field neither matches what has been seen in the 
evolution of the linearly polarized flux density during 
well-resolved, distinct events that have been followed in 
detail, nor is the spectral evolution of the position angle 
in the rise portion of the event generally consistent with 
the observations. Based on this evidence, it is concluded 
that the magnetic field within the emitting region cannot 
be largely of the ordered helical type. 

6. CONCLUSIONS 

The 'shock in jet' model for cm-waveband blazar vari- 
ability has been revisited, allowing for arbitrary shock 
orientation with respect to the jet fiow direction, and 
both random and ordered magnetic field. Oblique shocks 
can explain events with swings in polarization position 
angle much less than the 90° associated with transverse 
structures, while retaining the general characteristics of 
outbursts, including spectral behavior and level of peak 
percentage polarization. 

Specifically, as can be seen from Table [U the model 
reproduces the temporal behavior of the UMRAO data: 
the trend in spectral evolution through an opaque phase, 
while remaining quite fiat throughout the burst evolu- 
tion; significant levels of percentage polarization (but far 
less than the maximum permitted for a homogeneous 
source); the magnitude of temporal evolution seen in 
EVPA; and the frequency spread in EVPA during out- 
burst of order degrees. The range of models provided 
by varying the observer orientation with respect to the 
shock plane ( §5.2) Fig. [6|) display the same general char- 
acteristics, but a range of percentage polarization from 
almost zero to ^ 20%, corresponding to the range seen 
in UMRAO data (gUl]). 

For certain azimuthal orientations with respect to the 
shock normal the percentage polarization during out- 
burst is low - not much more than that seen in the quies- 
cent state - but in general, levels of polarization commen- 
surate with that seen in UMRAO monitoring data occur, 
indicating that oblique shocks can explain the data with- 
out recourse to special observer orientation. Similarly, 
for the fiow parameters adopted here (modest Lorentz 
factor, leading to modest aberration) oblique shocks can 
give rise to the observed levels of percentage polarization 
for a broad range of orientation with respect to the flow 
axis, as is the case for a transverse sho ck with simila r 
parameters. As found by Hughes, Aller fc AUeri (jl989b[ ). 
only modeling of the detailed outburst profile for spe- 
cific events, with the inclusion of retarded time effects 
(in addition to matching the precise level of polariza- 
tion achieved at peak outburst) will better constrain the 
viewing angle. 

Models dominated by a force-free, minimum energy 
magnetic field configuration (essentially helical) display 
a limited rise in percentage polarization, and a frequency 
dependent swing in polarization position angle not in 
agreement with the results of single-dish monitoring pro- 
grams, implying that the field cannot be dominated by 
an ordered component with the character of a force-free 



fiux rope in the quiescent state. This is in agreement 
with the conclusions of a number of studies (discussed 
in !j3.2|) that argue against a predominantly helical field 
configuration within the emitting region. 

Outbursts well-explained by the 'shock in jet' model 
are present during 7-ray fiaring in several sources, sup- 
porting the idea that shock events are associated with ac- 
tivity from the radio to 7-ray bands. Some 7-ray fiares 
might be associated not with propagating shocks, but 
rather with propagating particle density /magnetic field 
enhancements that enco unter stationary st ructures such 
as a recoUimation shock ()Agudo et al.l20l"ll) . Such events 
are best explained if the magnetic field of the propagat- 
ing enhancement is random, and so our conclusion that 
the quiescent flow of blazars is predominately random 
provides support for this explanation. 

We thank the referee for a careful reading of the 
paper, and many detailed, thoughtful, and construc- 
tive comments. This work was made possible by sup- 
port from NSF grant NSF0607523, NASA Fermi grants 
NNX09AU16G and NNX10AP16G, and by support for 
the operation of UMRAO from the University of Michi- 
gan. This research has made use of data from the MO- 
JAVE database th at is maintained by the MO JAVE team 
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Table 1 

Typical Flare Properties 



THse/T a P%max AEVPAt,„e AEVPA/^, 



A. UMRAO Data 



14.5 GHz 0.5 7 

8.0 0.6 4 

4.8 0.7 3.5 

30 - 110° 

start +0.30 20° 

max +0.5 

end -0.2 0° 



B. Shock Model A 

14.5 GHz 04 14 

8.0 0.4 12 

4.8 0.5 9 

25° 

start -0.11 4° 

max +0.13 1.5° 

end -0.13 4° 



C. Shock Model 01 (?? = 90°, 6*068 = 20°) 

14.5 GHz OA 14 

8.0 0.4 12.5 

4.8 0.4 9.5 

85° 

start -0.10 7° 

max +0.07 1.5° 

end -0.11 19° 



Note. — Specific values of a vary from event to event. Those listed 
show a typical change during an observed flare. 



Table 2 

Model Parameters 



Parameter 


Run A Runs B 


Runs C Runs D 


Semi-angle (fj.) 


2.4° 




Spectral index (a) 


0.25 




Fiducial LF (7^) 


1000.0 




Cutoff LF {7i) 


50.0 




Order multiple (/) 




0.1, 0.3, 1.0, 3.0 


Bulk LF (7y) 


2.5 




Sense 


F 




Shock width (w) 


0.4 




Compression (k) 


0.7 




Obliquity (rj) 


45° 


45, 90° 


Orientation (tp) 


0° 




Observer Ooi,^ 


10° 


20, 40, 60° 






30, 50, 70° 


Observer ipoba 


90° 0°, 45°, 90°, 135° 



180°, 225°, 270°, 315° 
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A = 4.8 GHz o = 8.0 GHz x = 14.5 GHz 



UMRAO Dec 4, 2010 
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Figure 1. Prom bottom to top: daily averages of the total flux density, the linearly polarized flux density, and the electric vector position 
angle in PKS 1510-089. The data at 14.5, 8.0 and 4.8 GHz are denoted by crosses, circles and triangles respectively. The radio band flaring 
shown is temporally associated with 7-ray activity. Note the differences in spectral behavior during the flaring periods shown. 
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Figure 2. From bottom to top: daily averages of the total flux density, the fractional linear polarization, and the electric vector position 
angle for OT 081 (1749+096). Symbols are as in Fig.[Tl 
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Figure 3. From bottom to top: 7 ray light curve and daily averages of the total flux density, the linearly polarized flux density, and 
the electric vector position angle for OJ 287 (0851+202). Symbols are as in Fig. [T] for the UMRAO monitoring data; source-integrated 
VLBA data obtained from the MOJAVE (15 GHz) and BU (43 GHz) websites are shown for comparison. The 7-ray light curve was kindly 
provided by S. Jorstad. Units are photons/sec/cm^xlO~ 
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A =4.8 GHz 



o =8.0 GHz 
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UMRAO Dec 6, 2010 
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Figure 4. From bottom to top: 60 day averages of the total flux density, the linear polarizati on, and the electric vector position angle for 
NRAO 530. Symbols are as in Fig. [T] A rotation measure of -63 rad/m^ has been assumed I IHuski n"988l). The Faraday-corrected EVPA 
is near 70° at 14.5 GHz during the relatively quiescent period since circa 2000 following a very large outburst in the 1990s. The position 
angle of the jet determined from MOJAVE measurements at the same frequency is 12° (jKovalev et al.ii2005i V 
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Figure 5. Oblique siiock evolution using the parameter set 'Run A' in Table [2] The panels and symbols correspond to those of Fig. [T] 
showing from bottom to top: the total flux density, the fractional linear polarization, and the electric vector position angle, at three 
frequencies corresponding to the UMRAO observations at 14.5, 8.0 and 4.8 GHz (crosses, circles and triangles respectively). 



17 



0.4 0.6 o.e 



0,4 0.6 O.e 

Time (years) 




0.4 0,6 O.e 

Time (years) 



0.3 0.4 0.6 

Time (years) 




0,3 0,4 0.6 0,8 1 



Time (years) 



0,6 o,e 



Time (years) 



Figure 6. The same model as shown in Fig. (5] for a range of observer orientations with respect to the shock plane. From panel (A) the 
azimuthal angle of the observer is 0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315°. (Note that in these models the zero point for the EVPA 
is arbitrary and can change from plot to plot — only the range in values is significant.) 
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Figure 7. Comparison of a transverse sliock model with that of oblique Run A, as shown in Fig. \5\ for a range of polar angles. In the 
left panels the shock orientation has been adjusted to be transverse to the flow, and in panels (A-C) the observer is viewing at angles 20°, 
40°, and 60° to the jet axis. The right panels are for the oblique case, and in panels (D-F) the observer is viewing at angles 30°, 50°, and 
70° to the jet axis. (The range of angles is chosen to optimally display variation with observer orientation. Note that in these models the 
zero point for the EVPA is arbitrary and can change from plot to plot - only the range in values is significant.) 
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Figure 8. The same model as shown in Fig. \5\ for a range of ordered to random magnetic field strengths. Prom panel (A) the 'order 
multiple' / is 0.1, 0.3, 1.0, and 3.0, so that in panel (A) the random field dominates, while in panel (D) the ordered field dominates. (Note 
that in these models the zero point for the EVPA is arbitrary and can change from plot to plot - only the range in values is significant.) 



